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Abstract 

We apply Stochastic Quantization Method to dissipative systems at finite tempera- 
ture. Especially, the relation of SQM to the Caldeira-Leggett model is clarified ensuring 
that the naive Wick rotation is improved in this context. We show that the Langevin 
system obtained by the "Improved Wick Rotation " prescription is equivalent to an ideal 
friction case ( low temperature limit) in the C-L model. We derive, based on our ap- 
proach, a general formula on the fluctuation-dissipation theorem for higher derivative 
frictions. 



^-mail address: nakazawa@nbivax.nbi.dk, nakazawa@psl.yukawa.kyoto-u.ac.jp 

2 Permanent address: Department of Physics, Shimane University, Matsue 690, Japan 



Quantization of dissipative systems has been discussed especially in connection with 
the description of non-equilibrium quantum phenomena. To describe the quantum effect 
in dissipative systems, one may naturally consider the dissipative system as a part of a 
"whole system ". In a physical point of view, the macroscopic dissipation comes from 
the microscopic interaction between the system and its environment. Realizing this idea 
explicitly, Cardeila and Leggett (C-L) developed a model to describe dissipation They 
introduced artificially harmonic oscillators as the environment, which couples with the 
system in question linearly. By integrating out the environmental degrees of freedom, 
they derived an effective action for the quantum effects with dissipation. 

In the previous paper 0, we applied Stochastic Quantization Method (SQM) || on 
the description of the quantum effect with dissipation. Since SQM can be formulated 
without action principle, namely it is formulated in terms of a Langevin equation which 
includes an equation of motion as the drift force, the artificial introduction of the envi- 
ronmental degrees of freedom is not necessary for the description of dissipative systems. 
This point may be an advantage of SQM in formulating a quantization of dissipative 
systems because it avoids the possible dependence of the physical results on the choice 
of the artificial environment and its interaction with the system in question. Our basic 
motivation is to clarify, whether SQM can describe a system interacting with its environ- 
ment which causes dissipation by starting from the phenomenological equation of motion 
including a dissipative term. In Ref. M, we considered the application of SQM to quan- 
tum mechanical systems with friction. We showed that some improvements are necessary 
to describe dissipative systems in SQM. As a guiding principle, we took our attention 
especially on the correspondence of our approach to the fluctuation-dissipation theorem 
and proposed an improved Langevin equation to describe the quantum mechanical sys- 
tems with dissipation at finite temperature. In this paper, we show the equivalence of the 
SQM approach in terms of the improved Langevin equation to the C-L model approach 
at finite temperature. 

We first recapitulate the essence of Ref. [2]. Let us consider the system with dissipation 
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for which the equation of motion is given by, 

d 

Mq + —V(q)+r]q = 0. (1) 

Here rj is a phenomenological friction coefficient. The dot denotes the derivative with 
respect to the real time coordinate t. In a naive sense, we may start with a Langevin 
equation^, 

^-q(r, t) = -i{Mq + ^-V(q) + V q} + £(r, t), (2) 
where the correlation of the white noise £ (r, t) is defined by 

< £(r, t)f (V, >= 2M(r - r')8{t - t'). (3) 

In the systems without dissipation, (2) with r\ = 0, the finite temperature Langevin equa- 
tion is derived in a heuristic way as follows. 

(a) : Continue the real ( Minkowski ) time to the imaginary ( Euclidean ) time by Wick 
rotation; t — ► —it. 

(b) : To improve the Langevin equation such that it has an equilibrium limit, "Wick " 
rotate the fictitious time r into — it. ( This procedure makes the drift force positive 
definite in the Langevin equation.) 

(c) : Redefine the white noise variable £ — > z£ to preserve the noise correlation (3) invari- 
ant under these continuations, (a) and (b). 

(d) : Require that all the variables satisfy the periodic ( anti-periodic ) boundary condi- 
tion for bosonic system ( fermionic system ). 

The procedure (a)-(d) really recovers the imaginary time method at finite temperature |Q || 
Following to this naive 'Wick rotation " prescription, we found that we couldn't obtain 



1 For the equation of motion (1), one can write a time dependent lagrangian. Then, we would start 
from the explicitly time dependent Langevin equation, 

JU(t, t) = ~if(t){Mq + ^-V(q) + w) + £(t, t), 

with f(t) = eO'/M")*. In this case, however, the explicit time dependence would cause some difficulties. 
It is also shown that the approach in terms of the above Langevin equation is, in a rigorous sense, not 
equivalent to the C-L approach)^]. 



the Langevin equation which describes the damped oscillator system at finite tempera- 
ture. As it is clarified in Ref. [2], the failure comes from the retarded ( causal ) nature 
of the Green's function obtained at zero temperature. Then we proposed the following 
improved Langevin equation for a damped oscillator case. 

^n(r) = -M{u 2 n + u)l)q n {r) - r]\u n \q n (r) + | n (r), 
< U{r)L(r') > = 2M(r - r')(3h5 m+nfl . 

(4) 

The expectation value, < >, means the finite temperature expectation value. The tilde 
denotes its Fourier component^. Here we have taken the additional procedure 
(e): Take the absolute value of the Fourier momentum in the friction dependent term in 
the momentum space Langevin equation. 

It was shown that the momentum space Langevin equation (4) really recovers the result 
which is consistent to the fluctuation-dissipation theorem. As it can be seen from (4), 
the absolute value prescription has been taken to make the drift force positive definite. 
We notice that the absolute value \uj n \ would not appear in the right hand side but u n 
would appear if we would consider the naive Wick rotation (a)-(d). Thus the absolute 
value prescription (e) for the friction dependent term is necessary for the existence of an 
equilibrium distribution at the infinite fictitious time limit. For a general potential, V(q), 
one can obtain the Langevin equation by simply replacing 

r f3h f) 

MulUr) - J o dte^—Viv + q(r, t)), (5) 

in the r.h.s. of (4). Here the vacuum v satisfies -§^V{v) = 0. 

Before showing the equivalence of our approach to the C-L model approach, in demon- 
strating how our approach works well, we consider the case of higher derivative frictions 
2 The coefficients of the Fourier series are defined by 

n 

with cu n — t denotes the imaginary time. 
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which is defined by the following phenomenological equation of motion. 

Mq + ^V(q) + (-)'- V 2 '^ = F ex , (6) 

where q^ 21 ^ = d 2l ~ 1 q(t)/dt 2l ~ 1 . The external force F ex {t) is introduced for the later 
convenience. / = 1 corresponds to Ohmic dissipation in (1) and I = 2 corresponds to 
super-Ohmic case|J. The sign factor, (— comes from that fact that the system is 
dissipative. Actually, the average of the energy, ^Mq 2 +V(q) , in a periodic motion is dis- 
sipated with the rate i]{d l q(t)/d l t} 2 . Thanks to the improved Wick rotation prescription 
(a)-(e), we may begin with the Langevin equation 

^Ur) = -M(u 2 n + co 2 )q n (r) - H^M 2 '" 1 ^) + | n (r), 

< £m(V)£n(V) > = 2%8{t - T)(3M m+nfi . 

(7) 

for F ex = 0. The factor (— comes from the same factor appears in (6). For simplicity, 
we mainly consider the harmonic oscillator with I = odd case, however, the generality is 
obviously recovered by the replacement (5) in (7). By using the solution of (7), we obtain 
the two-point function at finite temperature 

<q 2 > = lip <q(r,t)q(r',t)>, 

4r P { M(u; 2 + co 2 ) + (-y-i v \ Un \n-i /' 



(8) 



By using the same technique in Ref.[l], we obtain an equivalent expression to (8) as 
follows, 

h f°° . ,f3hu, 



2 n f°° . ,phuj, T . . , 

< q > = — coth( )lmx q {uJ)duj 

ix Jo 2 

E Im / ixM «2 A d", 



QpGu.h.p. 



M(u 2 -u 2 - 2i 1 u 21 - 1 ) ' 



7 = r]/2M, 



(9) 



where a p denotes the pole of Xq{ u ) i n the complex cu-plane. The summation J2 ap eu.h. P . 
means that the only poles on the upper half tu-plane should be summed up. For the case 
of I = 1, Ohmic dissipation, the formula turns out to be equivalent to the result appeared 
in Ref.[l]0 . We notice that the second term in the r.h.s. of the first relation in (9) 
is not a correction on the fluctuation-dissipation theorem which connects < q 2 > to the 
imaginary part of the response function , lmx q , but a precise reflection of the statement of 
this theorem that the contributions which come from only the retarded Green's function 
defined by the equation of motion (6) should be included in this formula. If the second 
term would not appear in the r.h.s. of (9), the causal structure of this formula would 
be violated. In this sense, the second term in (9) is recognized as the counter term for 
the contribution from u.h.p., possibly included in liax q , which may violate the causal 
structure of this formula. 

We remark that the Fokker-Plank equation which corresponds to (7) with (5) gives the 
following equilibrium distribution functional in momentum space at the infinite fictitious 
time limit, 

P cq [q] = iVexp[--]T — {Mul + R^KI 2 '- 1 }^! 2 - - jf dte^Viv + q(r,t)) , 

(10) 

provided that the potential ensures the positive definite drift force in (7) such as in the 
double well potential case||. For the harmonic oscillator with / = even, the Langevin 
equation (7) does not have any equilibrium limits. This is a difficulty not only for SQM 
approach but also for the corresponding C-L model approach because, as we will discuss 
later, the path-integral does not converge in this case. 

The remaining part of this paper is devoted to show the equivalence of the prescription 
(e) to the microscopic description in terms of C-L model. Now we clarify the relation of 

(7) with (5) to C-L model at finite temperature. The effective action of C-L model at 
3 In the Ohmic dissipation, 1=1, the second term is not necessary because the poles appears only on 
the lower half of the complex w-plane. 
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finite temperature is given by 

r Ph M rf3h rf3h 

S c s= dr[—q 2 + V{q)}+ dt dt'K{t -t')q{t)q{t'), (11) 
Jo I Jo Jo 

We notice that the dot denotes the derivative with respect to the imaginary time t in 
the following arguments. The Fourier components of the kernel function K(t) is defined 
byii 

K M = ^ ir^-^. (12) 



IX Jo UJUJ Z + UJ^ 

The explicit form of the spectral function J(uj), which is determined by requiring that 
the phenomenological equation of motion (1) should be recovered in a classical limit, is 
not crucial in the following arguments. This is because, from (12), the kernel function 
K(uj n ) is an even function of u n , a function of |a> n |, ensuring the validity of the absolute 
value prescription (e). For the Ohmic dissipation, it is required to be J{uS) = r\uj. While 
in a physically realistic situation, the dependence of J(uS) on the high frequency u would 
be modified depending on its microscopic physics. We here only assume, 

J{uj) = 7^e^ /LJc , (13) 

by introducing a critical value u c ^. Then K(uj n ) is given by 

KM = l j" n \, ■ (14) 
2 1 + \uj n \/uj c 

Now the correspondence of our approach and C-L model is obvious. At the limit u c — > oo, 
the Fokker-Planck distribution (10) is equivalent to path-integral measure via the effective 
action of C-L model. Namely, in this limit, (4) is equivalent to the Langevin equation 
derived from the effective action of C-L model at finite temperature, 

±q{T,t) = - S -^(T,t)+£{T,t). (15) 



4 In this effective action, the necessary subtraction of the "divergent " term to recover the phenomeno- 
logical equation of motion such as (1) and (6) is included]!]]^!?]]. Thus the potential appears in (11) is 
recognized as the renormalized one. The reader desiring a full account should consult Ref.p]. 
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Actually, we have the following Langevin equation as the configuration space represen- 
tation of (4)5 

>• « - " f - S P' ufSfe^ + <"> 

This correspondence is easily extended to the higher derivative friction cases. For 
the higher derivative case, to recover the phenomenological equation of motion (6), the 
spectral function is assumed to be 

J(w) = r^-V^, (17) 

Then we have the kernel function K(u n ) in the limit uj c — > oo, 

lim K(uo n ) = (-) l -^\uj n f-\ (18) 

This ensures the complete equivalence of the Fokker-Planck distribution (10) to the path- 
integral measure via C-L model effective action. 

In a physically interesting case, / = 2 ( super Ohmic case ), the explicit form of the 
Langevin equation in configuration space is given by 



^ ( t) M" ^ ^ ^ 
dr ' dq 7r dt 



^ dt ,_ q{r,t) - q{r,t') 



o 



+ f(r,t). (19) 



{(f)sin[vr(t-t')//%]} 2 . 
The friction dependent term in the r.h.s is the variation of the effective action of C-L 
model in super Ohmic case||. 

It is clear from (18) that the Langevin equation for the higher derivative friction case 
can not be obtained by the naive Wick rotation prescription from the friction dependent 
term rjd 2l ~ 1 q(t) / dt 21 ^ 1 . We have thus confirmed the validity of the absolute value pre- 
scription (e) in these well. We notice, however, the appearance of the factor 

(— ) in the phenomenological equation of motion (6), which comes from the positivity 
5 We also find another useful form of the configuration space expression of (4), 

—q(T,t)=Mq-- / dt — — h£(r,t). 



dr ' dq 7T J — oo (*-*') 

Here we define the integral by its principal value 



of the spectral function J(uS) as is clearly recognized in C-L model approach, would cause 
a difficulty not only in SQM approach but also in C-L model approach f\ In some cases 
with I = even, for example a harmonic oscillator, the path-integral does not converge 
with the effective action (11). The best way to derive an effective action in C-L model 
which avoids this difficulty would be to introduce a cut off in the high frequency part of 
the spectral function such as, 



where u c is a critical value of the frequency which is determined by the requirement for 
the effective action (11) to be bounded from the below. This would yield a modification 
on the phenomenological equation of motion such as (6) ensuring a positive definite drift 
force in the improved Langevin equation. It is, however, beyond C-L model approach as 
well as SQM approach to justify such an ad hoc modification (20). It depends on the 
detail of the microscopic physics and we do not discuss on this possibility further. In any 
cases, what we would like to remark is that a well-defined path-integral measure with the 
C-L model effective action which is bounded from the below corresponds to an improved 
Wick rotation prescription in SQM via the kernel, K(uj n ), the possible explicit form of 
which is fixed as a function of \u n \ by dimensional analysis up to a numerical factor. 

In this short note, we have clarified the equivalence of the "improved Wick rotation 
" prescription in SQM at finite temperature to C-L model approach. Our main claim is 
that, to derive the Langevin equation by starting from the phenomenological equation of 
motion, the naive Wick rotation should be improved. We also have derived a generalized 
formula corresponding to the fluctuation-dissipation theorem in higher derivative friction 
cases. We have mainly studied the limit that the low frequency dependence of the spectral 
function is extrapolated to the high frequency region {u c — > oo in (13) and (17)). In a 
physical point of view, this assumption corresponds to the low temperature case, i- << 
6 For the super Ohmic case, the opposite sign factor appears, as it is clear from (7), in the friction 
dependent term which comes from (6). In this case, if the friction coefficient rj would exceed its critical 
value which makes the friction term comparable to the "kinetic + potential " term, we could not obtain 
the equilibrium limit in the Langevin equation (7) with (5). 




(20) 



lo c . Thus it is more precise to conclude that the improved Wick rotation prescription in 
SQM in quantizing dissipative systems has its validity at low temperature. 
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